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Thermal conductivity in harmonic lattices with 
random collisions. 
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Olla 


Abstract We review recent rigorous mathematical results about the macroscopic be¬ 
haviour of harmonic chains with the dynamics perturbed by a random exchange of 
velocities between nearest neighbor particles. The random exchange models the ef¬ 
fects of nonlinearities of anharmonic chains and the resulting dynamics have similar 
macroscopic behaviour. In particular there is a superdiffusion of energy for unpinned 
acoustic chains. The corresponding evolution of the temperature profile is governed 
by a fractional heat equation. In non-acoustic chains we have normal diffusivity, 
even if momentum is conserved. 
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1 Introduction 

Lattice systems of coupled anharmonic oscillators have been widely used in order to 
understand the macroscopic transport of the energy, in particular the superdiffusive 
behavior in one and two dimensional unpinned chains. While a lot of numerical 
experiments and heuristic considerations have been made (cf. ll^[^l45l and many 
contributions in the present volume), very few mathematical rigourous scaling limits 
have been obtained until now. 

For harmonic chains it is possible to perform explicit computations, even in the 
stationary state driven by thermal boundaries (cf. ||43| ). But since these dynamics are 
completely integrable, the energy transport is purely ballistic and they do not pro¬ 
vide help in understanding the diffusive or superdiffusive behavior of anharmonic 
chains. 

The scattering effect of the non-linearities can be modeled by stochastic pertur¬ 
bations of the dynamics such that they conserve total momentum and total energy, 
like a random exchange of the velocities between the nearest neighbor particles. We 
will describe the results for the one-dimensional chains, and we will mention the 
results in the higher dimensions in Section |9] In particular we will prove how the 
transport through the fractional Laplacian, either asymmetric or symmetric, emerges 
from microscopic models. 

The infinite dynamics is described by the velocities and positions {{pxiqx) S 
of the particles. The formal Hamiltonian is given by 

^Y^pl + ]^Y^(^x-x^qxqx'. ( 1 ) 

X ^ x,x' 

where we assume that the masses are equal to 1 and that a is symmetric with a finite 
range or at most exponential decay |aj| < We define the Fourier transform 

of a function /: Z — >• R as /(k) = for k € T. We assume that a{k)>0 

for kfi^O. The function a){k) = xja{k) is called the dispersion relation of the chain. 


Unpinned Chains 

We are particularly interested in the unpinned chain, i.e. dc{0) = 0, when the total 
momentum is conserved even under the stochastic dynamics described below. Then, 
the infinite system is translation invariant under shift in q, and the correct coordi¬ 
nates are the interparticle distances (also called stretches, or strains); 

rx=qx-qx-\, xGZ. ( 2 ) 

When a"(0) > 0 we say that the chain is acoustic (i.e. there is a non-vanishing 
sound speed). We will see that this is a crucial condition for the superdiffusivity of 
the energy in one dimension. For unpinned acoustic chains we have that co{k) ~ \k\ 
as k ^0. 
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Pinned chains 

When a (0) > 0, the system is pinned and translation invariance is broken. In this 
case (o{k) ^ as k ^ 0. This is also the case for unpinned non-acoustic chains, and 
this is responsable for the diffusive behavior of the energy, cf. Section|7l 


Dynamics with Stochastic Collisions 

To the Hamiltonian dynamics we add random elastic collisions, where momenta of 
the nearest-neighbor particles are exchanged. This happens at independent random 
exponential times; each couple of particles labeled x,x+ 1 exchange their velocities 
Px and Pj^^i at exponential independent random times of intensity 7 . Equivalently 
there are independent Poisson processes {Nx^x+i {t), x G'Z} of intensity 7 , indepen¬ 
dent from the positions and velocities of all particles. The evolution of the system is 
described by the stochastic differential equations 


qx = Px 

Px = -{OC * q{t))x + {px+lit^) - Px{t^)) Nx,x+lit) + {px-lit^) - Pxit^)) Nx-l,x{t) 

(3) 

where Nx,x+i{f) = Tx.x+iU))’ with the random times when 

v +1 jumps, and px{t ) is the velocity of the particle x just before time f, i.e. 
lim 40 P;c(f-i). 

The evolution of the probability density on the configurations then follows the 
Fokker-Planck equations; 


dtf{t,p,q) = {A + S)f{t,p,q) (4) 

where A is the hamiltonian operator 

^ = E - i^qx’^)dpx) , (5) 

while S is the generator of the random exchanges 

Sf{p,q) = rE -/(p>?)) > ( 6 ) 

a: 

where is the configuration obtained exchanging px and Px+\- 

This stochastic perturbation of the Hamiltonian dynamics has the property to 
conserve the total energy, and in the unpinned case the resulting dynamics conserves 
also the total momentum {Y,xPx), and the volume or strain of the chain rx). It has 
also the property that these are the only conserved quantities. In this sense it gives 
the necessary ergodicity to the dynamics ( ll24l[T5]l I. 

We have also considered other type of conservative random dynamics, like a 
continuous random exchange of the momenta of each triplets {px-i,Px,Px+i}- The 
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intersection of the kinetic energy sphere j = C with the plane + 

Px + Px+\ = C gives a one dimensional circle. Then we define a dynamics on this 
circle by a standard Wiener process on the corresponding angle. This perturbation 
is locally more mixing, but it gives the same results for the macroscopic transport. 


Equilibrium stationary measures 

Due to the harmonicity of the interactions, the Gibbs equilibrium stationary mea¬ 
sures are Gaussians. Positions and momenta are independent and the distribution 
is parametrized, accordingto teh rules of statistical mechanics, by the temperature 
T = > 0. In the pinned case, they are formally given by 

Vp(dq,dp) - — - dqdp. 

In the unpinned case, the correct definition should involve the rx variables. Then 
the distribution of the r^’s is Gaussian and becomes uncorrelated in the case of the 
nearest neighbor interaction. For acoustic unpinned chains the Gibbs measures are 
parameterized by 


A = (j3 ' (temperature),p(velocity),T(tension)), 
and are given formally by 


Vx{dr4p) -- — - dr dp. (7) 

Non-acoustic chains are tensionless and the equilibrium measures have a different 
parametrization, see Section|2] 


Macroscopic space-time scales 

We will mostly concentrate on the acoustic unpinned case (except Section|7J. In this 
case the total hamiltonian can be written as 

d^{p^ ^) •= E y “ ^ L 

^ x,x' 

where qx — q^j = Lv=y+i ^ ^ ^define the energy per atom as: 

Px 1 / \2 

= y• 


ex{r,p) : 


(9) 
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There are three conserved (also called balanced) fields: the energy the 

momentum YxPx and the strain Yx Cc ■ We want to study the macroscopic evolution 
of the spatial distribution of these fields in a large space-time scale. We introduce 
a scale parameter e > 0 and, for any smooth test function 7 : R —R, define the 
empirical distribution 

e^7(ex)U;,(e-“f), = {r,,p,,ex). (10) 

.Y 

We are interested in the limit as e —0. The parameter a € [1,2] corresponds to 
different possible scalings. The value a = 1 corresponds to the hyperbolic scaling, 
while a = 2 corresponds to the diffusive scaling. The intermediate values 1 < a < 2 
pertain to the superdiffusive scales. 

The interest of the unpinned model is that there are three different macroscopic 
space-scales where we observe non-trivial behaviour of the chain: a = 1, |, 2. 


2 Hyperbolic scaling: the linear wave equation. 

Let us assume that the dynamics starts with a random initial distribution jig = (•)£ 
of finite energy of size i.e. for some positive constant Eq: 

£EM0))s < Eo. (11) 

where we denote ex(t) = ex(p(t),q(t)). We will also assume that some smooth 
macroscopic initial profiles 


uo(y) = (ro(y),Po(y),eo(y)) 


are associated with the initial distribution, in the sense that: 


lim^£ 

£^•0 


e£7(ex)U;,(0)- 



V5 >0, 


( 12 ) 


for any test function J. 

Then it can be proven (l28l) that these initial profiles are governed by the linear 
wave equation in the following sense 

= j J{y)'^{yd)dy (13) 

where u(y,f) = (r(y,f),p(y,f),e(y,f)) is the solution of: 

d,v=dyp, dtp = TidyV, a,e=Ti5y(pr) (14) 
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with Ti = -g^ (the square of the speed of sound of the chain). Notice that in the 
non-acoustic case a"(0) = 0, there is no evolution at the hyperbolic scale. 

Observe that the evolution of the fields of strain r and momentum p is au¬ 
tonomous of the energy field. Furthermore we can define the macroscopic mechan¬ 
ical energy as 

Cmech{y,t) = ^ (Tiv(y,tf-l-p(y,tf) (15) 

and the temperature profile or thermal energy as 

T{y,t) = e{y,t) - t,nech{y,t). (16) 

It follows immediately from (fT4l i that T(y,t) = T{y,Q), i.e. the temperature profile 
does not change on the hyperbolic space-time scale. 

There is a corresponding decomposition of the energy of the random initial con¬ 
figurations: long wavelengths (invisible for the exchange noise in the dynamics) 
contribute to the mechanical energy and they will evolve in this hyperbolic scale 
following the linear equations (fT4l i. This energy will eventually disperse to infin¬ 
ity at large time (in this scale). Because of the noise dynamics, short wavelength 
will contribute to the variance (temperature) of the distribution, and the correpond- 
ing profile does not evolve in this hyperbolic scale. See ll28l for the details of this 
decomposition. 


3 Superdiffusive evolution of the temperature profile. 

As we have seen in the previous section, for acoustic chains the mechanical part of 
the energy evolves ballistically in the hyperbolic scale and eventually it will disperse 
to infinity. Consequently when we look at the larger superdiffusive time scale £^"f, 
a > 1, we start only with the thermal profile of energy, while the strain and momen¬ 
tum profiles are equal to 0. It turns out that, for acoustic chains, the temperature 
profile evolves at the time scale correspondig to a = 3/2: 

limeY,J{£x) {exie^^^^t))£ = JJ{y)‘I{y,t)dy, t>0, (17) 

where ‘^{y^t) solves the fractional heat equation 

= -c\Ay\^l^% ^y,Q) = T{y,Q) (18) 

where c = This is proven in ED. We also have that the 

profiles for the other conserved quantities remain flat: 
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limeV 7(ex) (rv(e ^/^f))g = 0, 
e^o ^ 

\im eTj{ ex) {pjc{e^^^^t))e = Q, f > 0. 

c_i,n 


(19) 


Clearly the null value is due to the hnite energy assumption, otherwise it will be the 
corresponding constants, i.e. limg^o £Lx(Cr( 0 ))£ and lim£^o£L,v(/’j:( 0 ))£- 

In hnite volume with given boundary conditions (periodic or else), the mechani¬ 
cal energy will persist, oscillating in linear waves. At the larger superdiffusive time 
scale, waves will oscillate fast giving a weak convergence for the initial prohles of 
strain and momentum to constant values. 


4 The diffusive behavior of the phonon-modes. 

We have seen that ( fT9l l holds at the superdiffusive time scale and consequently at 
any larger time scale. But if we recenter the evolution of the strain and momentum 
around the propagation of the wave equation we see Gaussian Huctuations at the 
diffusive space time scale. For this purpose it is useful to introduce a microscopic 
approximation of the Riemann invariants (normal modes) of the wave equation: 



( 20 ) 


Once f^{t) are recentred on the Riemann invariants of the wave equation, they 
diffuse on the proper space-time scale, more precisely: 



For a proof see ll^ . 


5 Equilibrium fluctuations 

If we start with an equilibrium stationary measure corresponding to a certain tem¬ 
perature T, momentum p and strain r, then of course there will be no evolution 
of the empirical helds dehned by (fTOt . But we shall observe the evolution of the 
equilibrium time correlations dehned by 




(21) 
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where (•) denotes the expectation with respect to the dynamics in the corresponding 
equilibrium. Let us assume for simplicity of notation that r = 0 = p, otherwise we 
have to shift x along the characteristics of the linear wave equation. At time f = 0, it 
is easy to compute the limit (in a distributional sense); 

/TO ccT \ 

lime-'/25([e^V],0) = 5(y) 0 T 0 (22) 

\aT 0 {^ + a^)T^J 

where a is a constant depending only on the interaction. 

In the hyperbolic time scale a = 1 this correlation matrix evolves deterministi¬ 
cally, i.e. 

lime-i/2^([e-V],e^‘0=%,f) (23) 

£^0 

where 

= dyS'^^iyp), d,S^'^{yp) = (24) 

while for the energy correlations 

d,S^\y,t)=0. (25) 

In particular if p = f = 0, energy fluctuations do not evolve at the hyperbolic time 
scale. As for the energy profile out of equilibrium, the evolution is at a further time 
scale. By a duality argument (cf. EU), the evolution of the energy correlations oc¬ 
curs at the superdiffusive time scale with a = 3/2, namely 

= (26) 

e ->0 

where S^^{y,t) is the solution of; 


= -c\Ay\y^P\ 


(27) 


6 The phonon Boltzmann equation 

A way to understand the energy superdiffusion in the one-dimensional system is to 
analyze its kinetic limit, i.e. a limit for weak noise where the number of stochastic 
collisions per unit time remains bounded. In the non-linear case it corresponds to 
a weak non-linearity limit as proposed first in his seminal paper ll42l in 1929 by 
Peierls. He intended to compute thermal conductivity for insulators in analogy with 
the kinetic theory of gases. The main idea is that at low temperatures the lattice 
vibrations responsible for the energy transport can be described as a gas of inter¬ 
acting particles (phonons) characterized by a wave number k. The time-dependent 
distribution function of phonons solves a Boltzmann type equation. Over the last 
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years, starting from the work of Spohn Il44l . several papers are devoted to achieve 
phononic Boltzmann type equations from microscopic dynamics of oscillators. 

A rigorous derivation can be achieved for the chain of harmonic oscillators per¬ 
turbed by a stochastic exchange of velocities fT]. The main tool is the introduction 
of a Wigner function, which describes the energy density of the phonons. Let xjf be 
the complex field 


\(r{k) = ^{( 0 {k)q{k)+ip{k)), k G T, 

where p, q are the Fourier transform of the variables p, q and (o{k) is the dispersion 
relation of chain. The energy of the chain can be expressed in terms of the fields 
namely ^ The evolution of the field \j/ due to the the pure 

harmonic Hamiltonian without noise (7=0) reads 

d,^{k,t) = -ico{k)T^{k,t), 

therefore the quantities | vr(k)p are preserved by the harmonic dynamics. 

The Wigner distribution is defined in analogy to the usual one in quantum me¬ 
chanics 


W^iyXt) = (e/2) / 

m (28) 

W^^,k,t) := {xif{k-e^/2,t)* vf(k + e^/2,f)), 

where (•)£ denotes the expectation value with respect to the initial measure, chosen 
in such a way that the average of the total energy,is of order e^*, i.e. e {J^)e < Eq, 
see dnii- We also require that the averages of all px and qx are zero. 

The Wigner distribution W^{y,k,t) defined by (|28]) gives a different energy dis¬ 
tribution from the one considered in the previous sections, i.e. < ejg iyj (f) >. But in 
the macroscopic limit, as e —0 , they are equivalent lIZTll . 

We look at the evolution of the Wigner function on a time scale e^^t, with the 
strength of the noise of order e, i.e. we consider the dynamics defined by Q with 
/replaced by e/, like in a Boltzmann-Grad limit. This evolution is not autonomous 
and is given by 

drW^{^,k,t)=-i^ (O'{k)W^{^,k,t) + 7<^W" {^,k,t) 

- I, k,f) -1 ^,k,f)* + 

where is the field ,k,f) = {\(r(k + e^/2,t)Y{k — e^/2,t))e and*^ is a linear 
operator. The transport term is due to the harmonic hamiltonian, while the ’’colli¬ 
sion” operator ^ is related to the stochastic noise. 

It turns out that the field is the Wigner distribution associated to the difference 
between the kinetic energy and the potential energy. This is fast oscillating on the 
time scale and in the limit e —> 0 it disappears after time integration. Therefore 
in the limit e —> 0 the Wigner function W® weakly converges to the solution of the 
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following linear Boltzmann equation 

dtW{y,k,t) + (k) d^W{y,k,t) =^W(y,k,t). (30) 

2n 

It describes the evolution of the energy density distribution, over the physical space 
R, of the phonons, characterized by a wave number k and traveling with velocity 
Co'{k). We remark that for the unpinned acoustic chains co'{k) remains strictly posi¬ 
tive for small k. The collision term has the following expression 

^m= I dk'R{k,k'mk')-f{k)], 

Jr 

where the kernel R is positive and symmetric. One can write an exact expression on 
R, nevertheless its crucial feature is that R behaves like k^ for small k, due to the 
fact that the noise preserves the total momentum. Naively, it means that phonons 
with small wave numbers travel with a finite velocity, but they have low probabil¬ 
ity to be scattered, thus their mean free paths have a macroscopic length (ballistic 
transport). This intuitive picture has an exact statement in the probabilistic interpre¬ 
tation of (l30l l. The equation describes the evolution of the probability density of a 
Markov process {K{t),Y{t)) on T x R, where K{t) is a reversible jump process and 
Y (f) is an additive functional of K, namely Y (f) = Jq (o'{Ks)ds. A phonon with the 
wave number k waits in its state an exponentially distributed random time T(k) with 
mean value ~ k^^ for small k. Then it jumps to another state k' with probability 
~ k'^dk'. The additive functional Y (f) describes the position of the phonon and can 
be expressed as 

T(f) = f 

i=i 

where is the Markov chain given by the sequence of the states visited by 

the process K{t). Here denotes the number of jumps up to the time t. The tail 
distribution of the random variables {T(X,)a)'(X, )} with respect to the stationary 
measure n of the chain behaves like 

7t[\T{Xi)co'{Xi)\> X] ~ 

Therefore the variables T{Xi)(o'{Xi) have an infinite variance with respect to the 
stationary measure. We remark that the variance is exactly the expression of the 
thermal conductivity obtained in |[3l. The rescaled process N^'^/^YiNt) converges 
in distribution to a stable symmetric Levy process with index 3/2 ( ll26ll . H). As 
a corollary, the rescaled solution of the Boltzmann equation W{N^^^y,k,Nt) con¬ 
verges, as A —?► - 1 - 00 , to the solution of the fractional diffusion equation 

dtu{y,t) = -\A\^^‘^u{y,t). 

A different, more analytic approach can be found in 1^ . 
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In the pinned or the non-acoustic cases, co' ^ k for small k, T{Xi)co'{Xi) has finite 
variance with respect to the stationary measure. In particular this variance coincides 
with the thermal diffusivity computed by the Green-Kubo formula as. Then one 
can prove that the rescaled solution W{N^/'^y,k,Nt) converges to the heat equation 

d,u{y,t) =DAu{y,t), 
with D given by the thermal diffusivity. 

The results described above give a two step approach to the diffusion or superdif¬ 
fusion of the energy: first a kinetic limit where the Boltzmann phonon equation is 
obtained in the weak noise limit, then a superdiffusive or diffusive rescaling of the 
solution of this equation. The results described in Section| 6 ]concern a simple space- 
time rescaling, without any weak noise approximatiorQ. Still Boltzmann equation 
helps to understand the proof of (fTTl i. that goes under the following lines. 

Let us just consider the superdiffusive case and consider the evolution of the 
Wigner distribution at time and for the noise of intensity 7 . Equation ( |29] | 

becomes 


d,W‘^{^,k,t)=-i^(o'{k)e -fye 

70 1 \ 

This looks like a very singular limit. Still due to fast oscillations the T® terms 
disappears after time integration. Furthermore, because the number of collisions 
per unit time tends to infinity, the limit of the Wigner distribution homogenizes in 
the variable k, i.e. its limit becomes a function W(y,t) independent of variable k. 
Assume that the above facts have been proven, and consider the case of the simple 
random exchange of the velocities, that give a scattering rate of the form: R{k,k') = 
R{k)R{k'), with / R{k)dk = 1. The argument below, that follows the line of 1^ . can 
be generalized to various rate functions Gzl. The Laplace transform in time of 
the Wigner distribution satisfies the equation: 


(^X+lye ^^^R{k) + ico'{k)^e W£(A,^,k) 

= We{^,k,0) + e-^^^2YR{k) JR{k')dk'+ 


(32) 


By dividing the above expression by Dg{X,^,k) = e^/'^X+2YR{k) + ico' [k)^e, mul¬ 
tiplying by 2YR{k) and integrating in k one obtains 


ae j We{X,^,k)R{k)dk- J 


2YR{k)Ws{^,k,0) 

DeiX,^,k) 


dk = ff{e) 


(33) 


* In the non-linear cases we cannot expect that the two step approach would give the same result 
of the direct rescaling of the dynamics. In the j8-FPU the kinetic limit seems to give a different 
superdiffusion scaling than the direct limit I38II45I . 
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where, since R{k) ^ and due to the assumptions made on co'{k): 


= lye 


-3/2 


(-/ 


2yR^{k) 


dk\ —1- X-\-c\^ 


Thanks to the homogenization property of We 

f We{X,^,k)R{k)dk —1 w(A,^). 

J £^-0 

Furthermore 

2yR{k) J 

e^O 

and we conclude that the limit function w{X,^) satisfies the equation 

(A+c|^p/2)w(A,^) = jw{^X0)dk (34) 

where W(^,k,0) is the limit of the initial condition. Equation ( l34b is the Laplace- 
Fourier transform of the fractional heat equation. Assuming instead that (£)'{k) ^ k, 
a similar argument gives the normal heat equation lIZTll . 


7 Non-acoustic chains: beam dynamics 


We have seen in the previous sections that the ballistic behaviour at the hyperbolic 
scale and the superdiffusive behavior of the energy strictly depends on the positivity 
of the sound velocity Ti. In the case a{0) = a"{0) = 0 all the coefficients in these 
evolutions are null and in fact the limit of the energy follows a regular diffusion. 
Notice that the dynamics is still momentum conserving. 

A typical example is given by the following choice of the interaction: 


00 = 3/2, Oi=-l, 02 = 1/4, Ox = 0, x>2 (35) 

(then a{k) = 4sin'^(7rk)) that corresponds to the Hamiltonian 


^=E 




y + g i^x+i — '^qx + qx-if 
(36)' 


Notice that the expression (|2l) is not defined for any value of the tension T, this is 
why we also call these chains tensionless. Basically, the equilibrium energy does not 
change by pulling the chain. It does change by bending it, this is why the relevant 
quantities are defined by the local curvatures or deflections: 


Kx = -Aqx = 2qx - qx+l - qx-\- 


(37) 
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The relevant balanced quantities are now {Kx,Px,ex)- 

The invariant equilibrium measures are formally given by 


g-li[J^-p'LxPx-S,I.x 'Ta-I 


dtcdp 


(38) 


where the parameter £ is called load. Notice that these measures would be non¬ 
translation invariant in the coordinates Tv’s. 

Under these conditions the sound velocity Ti is always null, and there is no bal¬ 
listic evolution of the chain. In fact it turns out that the macroscopic evolution of the 
three quantities {i{t,y),p{t,y),c{t,y)) is diffusive. By defining 



after the corresponding space-time scaling we obtain ( 12911 ): 


(d,t=- d^p, 

, d,p = \d^t + Yd^p, 

5fT = Dyd^^I+l {dypf-pd^p 


(39) 


The first two equations are the damped Euler-Bernoulli beam equations. The third 
one describes the diffusive behavior of the energy. The thermal diffusivity Dy can be 
computed explicitely and diverges as 7—0 (the deterministic dynamics has ballistic 
energy transport as every harmonic chain), (the damping terms involving 7 are due 
to the exchange noise). In particular for constant initial values of f and p, the energy 
(temperature) profile follows a normal heat equation with the thermal diffusivity 
that can be computed explicitly (cf. l29l ). 

These models provide rigorous counter-examples to the usual conjecture that 
the momentum conservation in one dimension always implies superdiffusivity of 
the energy (cf l34l 1^ ). The presence of a non-vanishing sound velocity seems a 
necessary condition. 


8 A simpler model with two conserved quantities 

In this section we consider the nearest neighbors unpinned harmonic chain (with 
mass I and coupling forces ai = a_i = — ■^) but we add different stochastic colli¬ 
sions with the properties that they conserve the total energy and some extra quantity 
(that we call ’’volume”) but no longer momentum and stretch. By defining a = 
and the field { tJj: £ M ; x G Z} by 7]2x = opx and = px, the Hamiltonian equa¬ 
tions are reduced to 


rix = a{rix+i-rix-i). 
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The stochastic collisions are such that at random times given by independent Pois¬ 
son clocks Nx,x +\(0 of intensity 7 the kinetic energy at site x is exchanged with the 
corresponding potential energy. The simplest way to do it is to exchange the vari- 

n2 

able r\x with rix+ 1 . Because of the form of the noise the total energy ^ and the 
’’volume” Y.x'Hx = Lx(Pv + ari) are the only conserved quantities of the dynamics 
(El). By reducing the number of conserved quantities from 3 (energy, momentum, 
stretch) to 2 (energy, volume) we expect to see easily the influence of the other con¬ 
served quantity on the superdiffusion of energy. The nature of the superdiffusion 
for models with two conserved quantities are studied in the nonlinear fluctuating 
hydrodynamics framework by Spohn and Stoltz in BtI . 

In the hyperbolic time scaling, starting from an initial distribution associated 
to some smooth macroscopic initial volume-energy profiles (^ 0 ( 7 ), 60 ( 7 )), we can 
prove that these initial profiles evolve following the linear wave equation (v'(f ,7), e(f ,7)) 
which is solution of 

diV = 2adyV, dtZ = ady{v^). 

As in Section|2]we can introduce the mechanical energy Cmechiy,t) = 2 '’'^ the 

thermal energy T(7, f) = e( 7 ,f) ~ ^mcc/i(7T). The later remains constant in time. 

Mutatis mutandis the discussion of Sections [3 H] and |3 can be applied to this 
model with two conserved quantities with very similar conclusions El (MED- 
The interesting difference is that in (fTsT l and (l27T i the fractional Laplacian has to be 
replaced by the skew fractional Laplacian: 

d,‘l=-c{\Ay\^l^-Vy\Ay\^l^}%, ^(y, 0 ) = T (7, 0) (40) 

for a suitable explicit constant c > 0. The skewness is produced here by the inter¬ 
action of the (unique) sound mode with the heat mode. In the models of Section [T] 
which conserve three quantities, there are two sound modes with opposite velocities. 
The skewness produced by each of them is exactly counterbalanced by the other one 
so that it is not seen in the final equations (fTSl l and (l27l) . 


The extension problem for the skew-fractional Laplacian 

For the model with two conserved quantities introduced in this section, the deriva¬ 
tion of the skew fractional heat equation, at least at the level of the fluctuations in 
equilibrium as defined in Section |5] can be implemented by means of the so-called 
extension problem for the fractional Laplacian ii,El. As we will see, this exten¬ 
sion problem does not only provides a different derivation, but it also clarifies the 
role of the other (fast) conservation law (i.e. the volume). It can be checked that for 
any j3 > 0 and any p e R, the product measure with Gaussian marginals of mean p 
and variance (temperature) are stationary under the dynamics of {rix(t);x G Z}. 

Let us assume that the dynamics starts from a stationary state. For simplicity we as- 


2 In El only the equilibrium fluctuations are considered but the methods developed in 1271 . 1281 
can be applied also to the models considered in this section. 




Conductivity of harmonic lattices 


15 


sume p = 0 and a=\. The space-time energy correlation function f) is defined 

here by 

It turns out that the energy fluctuations are driven by volume correlations. Therefore 
it makes sense to define the volume correlation function as 

Ge{x,y,t) = {r]^{e-^/h)r]y{e-^/h) {rio{of -. 

Let / : [0,7’] xK.—^-Rbe a smooth, regular function and for each t G [0,7’], let 
M, : M X R+ —>■ M be the solution of the boundary-value problem 

r —dxU + ydyU = 0 
\ dyu{x,0) = dxf{x,t). 

It turns out that dxUf{x,Q) = ^ft{x), where ^ is the skew fractional Laplacian 
defined in ( l40l ). 

For test functions / : R —>■ R and m : R x R+ R define 
{Se{t)J)E = Y,SE{x,t)f{ex), 

jfGZ 

{Ge{t),u)E= Ge{x,y,t)u{^(x + y),y/e(j-x)). 

x,yel, 

x<y 

After an explicit calculation we have; 

)/?)e ~ {GE{t),Ut)E + £{SE{t),Ut{-,0))£'^ = {Sg{t),dxft)E 

plus error terms that vanish as e 0. From this observation it is not very difficult 
to obtain that, for any smooth function / of compact support, 

lim(5£(f),/)£ = j P{t,x)f{x)dx, 

where P{t,x) is the fundamental solution of the skew fractional heat equation (l40l i. 

Let us explain in more details why the introduction of the test function 
solves the equation for the energy correlation function Se{x,t). It is reasonable to 
parametrize volume correlations by its distance to the diagonal x = y. The micro¬ 
scopic current associated to the energy rix{t)^ is equal to rjx{t)rix+i{t), which can be 
understood as the volume correlations around the diagonal x =y. Volume evolves 
in the hyperbolic scale with speed 2. Fluctuations around this transport evolution 
appear in the diffusive scale and are governed by a diffusion equation. This means 
that at the hyperbolic scale £^', fluctuations are of order explaining the non¬ 

isotropic space scaling introduced in the definition of {GE{t),u)g. It turns out that 
the couple {Sg,Gg) satisfies a closed system of equations, which can be checked to 
be a semidiscrete approximation of the system 
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{ ^/ed,ut = -d^ut + ydyUt 
dyUt{x,Q) = dj,{x) 

dtft = d^Ut{x,0). 

Therefore, the volume serves as a fast variable for the evolution of the energy, which 
corresponds to a slow variable. The extension problem plays the role of the cell 
problem for the homogenization of this fast-slow system of evolutions. 


9 The dynamics in higher dimension 

One of the interesting features of the harmonic dynamics with energy and momen¬ 
tum conservative noise is that they reproduce, at least qualitatively, the expected 
behavior of the non-linear dynamics, also in higher dimensions. In particular in the 
three or higher dimension the thermal conductivity, computed by the Green-Kubo 
formula, is finite, while it diverges logarithmically in two dimension (always for 
non-acoustic systems), cf. EE). 

In dimension c/ > 3 it can be also proven that equilibrium fluctuations evolve 
diffusively, i.e. that the asymptotic correlation S^^{y,t), defined in section|7]but with 
a diffusive scaling, satisfies 0; 


d,S^^=Dd^P^ (41) 

where D can be computed explicitely by Green-Kubo formula in terms of CO and 
scattering rate R Em. Similar finite diffusivity and diffusive evolution of the fluc¬ 
tuations are proven for pinned models (6:(0) > 0), see 0. In two dimensions, while 
the logarithmic divergence of the Green-Kubo expression of the thermal conductiv¬ 
ity is proven in 0, the corresponding diffusive behaviour at the logarithmic cor¬ 
rected time scale is still an open problem. For the result obtained from the kinetic 
equation see 0. The two-dimensional model is particularly interesting in light of 
the large thermal conductivity measured experimentally on graphene ifSTI . an essen¬ 
tially a two-dimensional material. 


10 Thermal boundary conditions and the non-equilibrium 
stationary states 

Traditionally the problem of thermal conductivity has been approached by consid¬ 
ering the stationary non-equilibrium state for a finite system in contact with heat 
bath at different temperatures (cf. ll43l l34l ). This set-up is particularly suitable for 
numerical simulations and convenient because it gives an straight operational defi¬ 
nition of the thermal conductivity in terms of the stationary flux of energy, avoiding 
to specify the macroscopic evolution equation. But for the theoretical understanding 
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and corresponding mathematical proofs of the thermal conductivity phenomena, this 
is much harder than the non-stationary approach described in the previous section. 
This because the stationary state conceals the space-time scale. 

The finite system consists of 2N + 1 atoms, labelled by x = —N^ ■■■ ,N, with end 
points connected to two heat baths at temperature 7) and T^. These baths are modeled 
by Langevin stochastic dynamics, so that the evolution equations are given by 

qx=Px, x = 

px = -{Un * q{t))x + {px+\{t^) - Px{t^)) Nx,x+\{t) + {px-\{t^) - Px{t^)) Nx-\,x{t) 
+ ^x-N (y—P-N + \/^l dW-N{t + dx^N (^—PN + dWN{t , 

(42) 


where w-N(t),WN{t) two independent standard Brownian motions, and the cou¬ 
pling Un is properly defined in order to take into account the boundary conditions. 
For this finite dynamics there is a (non-equilibrium) unique stationary state, where 
the energy flow from the hot to the cold side. Observe that because of the exchange 
noise between the atoms, the stationary state is not Gaussian, unlike in the case 
studied in ll43l . 

Denoting the stationary energy flux by J^, the thermal conductivity of the finite 
chain is defined as 


Kjv 


lim 

|7i-TK0 


{2N+\)Jn 

Tl-Tr 


(43) 


For a finite N it is not hard to prove that Km can be expressed in terms of the cor¬ 
responding Green-Kubo formula. For the periodic boundary unpinned acoustic case 
this identification gives Kjv (cf. ||2l). For the noise that conserves only the 

energy, but not momentum (like independent random flips of the signs of the mo¬ 
menta), the system has a finite thermal diffusivity and the limit Kn ^ K,asN ^ + 0 °, 
can be computed explicitely, as proven in HI. 

The natural question is about the macroscopic evolution of the temperature pro¬ 
file in a non-stationary situation, and the corresponding stationary profile. It turns 
out that this macroscopic equation is given by a fractional heat equation similar to 
Ct]) with a proper definition of the fractional laplacian on the interval [—1,1] 

subject to the boundary conditions ‘T(—1) = 7;,T(1) = T,-. This is defined by using 
the following orthonormal basis of functions on the interval [—1,1]: 


M„(y) =cos(«a:(y-f l)/2), n = 0,1,... (44) 


Any continuos function f{y) on [—1,1] can be expressed in terms of a series ex¬ 
pansion in u„. Then we define \A\'^u„{y) = {nn/2)'^^Un{y). Observe that for i = 1 
we recover the usual definition of the Laplacian. For s ^ I this is not equivalent to 
other definitions of the fractional Laplacian in a bounded interval, e.g. ll5^ . 

Correspondingly the stationary profile is given by 
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|4rT = 0 , in (- 1 , 1 ) 
T(-1) = 7),5:(1) = 7;, 


(45) 


namely T(y) = \{Ti + Tr) + j{Ti- Tr)9{y), with 



m odd 


Cs such that 0(—1) = 1- This expression corresponds with the formula computed 
directly in ||3^ using a continuous approximation of the covariance matrix of the 
stationary state for the dynamics with fixed boundaries. 


11 The non-linear chain 

From the above rigorous results on the harmonic chain with the random collision 
dynamics, and the arguments of Spohn from fluctuating hydrodynamics and mode 
couplings (cf. Il45l and Hhl l. we can conjecture the corresponding behaviour in the 
anharmonic case. 

Let us consider just nearest neighbor interaction given by the potential energy 
V{qx — qx-i) of an anharmonic spring. We assume V ; R —?► (0, + 0 °) is smooth and 
that it grows quadratically at infinity. Define the energy of the oscillator x as 



(46) 


The dynamics is defined as the solution of the Newton equations 


qx{t)=px, Px{t) = -{V'{rx)-V'{rx+i)), xSZ 


(47) 


plus a random exchange of velocities as in the previous sections, regulated by an 
intensity 7 . The equilibrium Gibbs measures are parametrized by 


1 = (j 3 ' (temperature),/? (velocity ),t (tension)) 


and are given explicitly by 



(48) 


X 


When a random exchange of velocity is present (7 > 0) it can be proven that these 
are the only regular translation invariant stationary measures ( ll24lfT5]l l. We have 
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These thermodynamical relations can be inverted to express the parameters (j3^ *, p, T 
in terms of (p, r, e). It turns out that the tension is then equal to a nonlinear function 
f{r,u) of the average stretch f and the average internal energy S = e — p^/2. 

After the hyperbolic rescaling of the dynamics, we expect that the empirical dis¬ 
tribution of the balanced quantities converge to the system of hyperbolic equations: 

dtr = dyp 

dtp = dyf (49) 

dte = dy (pf). 

This limit can be proven, under certain condition on the boundaries, in the smooth 
regime, if the microscopic dynamics is perturbed by a random exchange of velocities 
between nearest neighbors particles particles by using relative entropy methods ( BTl 
1^1151 ). The limit should be still valid after shocks develop, with the limit profile 
given by an entropic weak solution. This is a main open problem in hydrodynamic 
limits. 

After a long time, the (entropy) solution of (l49l l should converge (maybe in a 
weak sense) to some mechanical equilibrium characterized by: 

p(y) =po, f (f(y),M(y)) = To, for some constants po, Tq. (50) 

It is very hard to characterize all possible stationary solutions that satisfy (l50l) . Prob¬ 
ably they are generically very irregular. But certainly if we start with a smooth initial 
condition that satisfies (fSOl l. they do not move. Also by the relative entropy methods, 
it is possible to prove that starting with such initial profiles, the empirical distribu¬ 
tion of the balanced quantities will converge at the hyperbolic space-time scale to 
such a stationary solution at any time. 

Still we do know that the microscopic dynamics will converge to a global equi¬ 
librium, so this implies that there exists a larger time scale such that these profiles 
will evolve and eventually reach also thermal equilibrium. 

There is a numerical evidence and heuristic arguments about the divergence of 
the Green-Kubo formula defining the thermal diffusivity for such one dimensional 
systems, so we expect that the larger time scale at which these profiles evolve is 
superdiffusive. 

From the nonlinear fluctuation hydrodynamics ( ll45ll ). one can conjecture the 
following: the space-time scale is (e^'x, and the temperature T{x,t) = 

evolves following some fractional heat equation, possibly non-linear. If 
V is symmetric and Tq = 0, then a = 3/4, and in all other cases a = 5/6. 


12 The disordered chain 

The effect of disorder on transport and phonons localization properties in chains of 
oscillators has attracted a lot of interest (||20l [HJ). Randomness can appear at the 
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level of the masses of particles or at the level of the potentials. We consider only 
the case of random masses with non random potential V or the case of non random 
masses and non random interaction potential V with random harmonic pinnings. In 
the first case, the Hamiltonian is then given by 

X ^ X,x' 

where {nix} are positive random variables, while in the second case, the Hamiltonian 
is given by 

xy X 

where {Vx} are positive random variables, m > 0 being the mass of the particles. 
The presence of randomness is relevant for the thermal properties of the system but 
the fact that randomness affect potentials or masses is not. 

For one dimensional unpinned disordered harmonic chains it is known that the 
behavior of the conductivity is very sensitive to the boundary conditions since it 
can diverge as ^/N or vanish as 1 /^/N with the systems length N (lfn [T9ll49ll50l l. 
If harmonic pinning is added localization of normal modes leads an exponential 
decay of the heat current and a zero conductivity. The situation in higher dimensions, 
even in the case of harmonic interactions, is still under debate but it is expected 
that conductivity is finite in dimension c/ > 3 if disorder is sufficiently weak (1321). 
About the effect of nonlinearities, numerical evidences suggest that a very small 
amount of anharmonicity in pinned chains is sufficient to restore a diffusive regime 
with a positive finite value of the conductivity ( EnH . However it is a challenging 
open question to decide if the transition from an insulator to a conductor occurs at 
zero or some finite small value of anharmonicity (ll9l 12711401 1. 

It is suggestive to think that a stochastic noise could affect transport properties of 
harmonic chains in some rough sense similar to the addition of nonlinearities. This 
question has been first address in nni in the Green-Kubo formula setting, revisited 
in II 22 I from the non equilibrium stationary state point of view (see section ITOli and 
extended in ifTTl l25l to incorporate weakly nonlinear chains. In these papers, the 
authors consider a disordered harmonic chain, or weakly nonlinear in lfT3ll25l . with 
a stochastic noise which consists to flip, independently for each particle, at indepen¬ 
dent random exponential times of mean 1/A, A >0, the velocity of the particle. No¬ 
tice that this energy conserving noise is very different from the noise considered in 
the rest of the paper since it does not conserve momentum. In particular, for ordered 
pinned and unpinned nonlinear chains, this noise is sufficient to provide a finite con¬ 
ductivity K ca. However, it turns out that for an ordered harmonic chain, ff ~ A ^ 
and, as expected, increases to infinity as the strength of the noise A / 0. In ifTTl l25l 
it is proved that localization effects persist: K= ^{X) for a pinned disordered chain 
with a small anharmonic potential, and tc ~ A for a pinned harmonic chain. As far 
as we know disordered chains with energy-momentum conserving noise have never 
been investigated. 
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